Impact of relativistic chiral one-pion exchange on nuclear matter properties 
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We study nuclear matter properties in a model that features relativistic chiral one-pion exchange 
and contact interactions between nucleons. We apply thermal field theory methods to compute 
microscopically the nucleon self-energy arising from one-pion exchange and its contribution to the 
energy per nucleon in isospin-symmetric nuclear matter as well as neutron matter. The contact 
interactions are fitted to reproduce the energy density of isospin-symmetric nuclear matter at sat- 
uration as well as benchmark calculations for neutron matter. We find that relativistic one-pion 
exchange contributes about half of the binding energy of nuclear matter, the remainder being pro- 
vided by the contact terms. We compare our model results to non-relativistic calculations based on 
low-momentum nucleon-nucleon potentials and discuss the symmetry energy and the compressibility 
of isospin-symmetric nuclear and neutron matter. 



I. INTRODUCTION 



The equation of state of bulk nuclear matter has at- 
tracted considerable attention over time, as it has a sub- 
stantial impact on the properties of neutron stars as well 
as finite nuclei. The equation of state is determined 
by two, generally related, ingredients: the force be- 
tween nucleons and the many-body approximation used 
to compute the thermodynamic properties of nuclear 
matter. Nuclear interactions are accurately modeled in 
terms of potentials arising from meson exchange: on the 
one hand models have been used which are based on 
phase-shift equivalent one-boson exchange [l[ and, on 
the other hand, more recently, (multi)-pion-exchange po- 
tentials based on chiral power counting Q. The diver- 
sity of ab-initio many-body approaches includes varia- 
tional Monte Carlo methods [3[ and propagator meth- 
ods with (medium-) rcnormalized soft interactions as, 
e.g., in non-relativistic and covariant Brueckncr-type the- 
ories Hartree-Fock calculations of nuclear matter 
properties were revived with the advent of soft chiral po- 
tentials Q which, however, require three- and four-body 
forces to account for contributions already included at 
the two-body level in the non-perturbative resummation 
schemes. 

In this work we apply thermal field theory (TFT) 
methods to compute the equation of state of isospin- 
symmetric nuclear and neutron matter. The method has 
been applied in the past in the context of QED and QCD 
plasmas to compute the quasiparticle energies of elec- 
trons and quarks, respectively, in a thermal medium Q. 
In the low-energy regime of interest to us, the relevant 
degrees of freedom are nucleons and pions. Chiral sym- 
metry is an (approximate) symmetry of the strong in- 
teraction. Because the chiral symmetry is spontaneously 
broken in Nature, pions emerge as the (pseudo-) Gold- 
stone bosons of the theory. The description of the strong 
interaction in terms of chiral Lagrangians admits cer- 
tain approximation schemes in terms of power counting 
of small quantities 0, Q > which in principle allow for a 



systematic order-by-order improvement of a given calcu- 
lation. Below, we will combine relativistic TFT methods 
and the chiral Lagrangian description of nuclear inter- 
actions to address the computation of thermodynamic 
properties of nuclear and neutron matter. 

Approaches similar to ours were previously developed 
in Refs. [§-[l2|. Fraga et al. Q derived analytical expres- 
sions for the zero-temperature self-energy of a nucleon 
due to pion exchange in dilute nuclear matter to leading 
order in the chiral expansion (but they did not address 
nuclear saturation). Lutz et al. and Kaiser et al. [ll| 
used chiral Lagrangians and expansions in small Fermi 
momentum to construct equations of state in the heavy- 
baryon limit. The saturation in Ref. fiol ] arises due to 
correlations induced by the one-pion-exchange interac- 
tion. In Ref. 11( two-pion exchange produces nuclear 
binding at three-loop level with a suitably adjusted mo- 
mentum cut-off. Effective field theory of nuclear matter 
with nucleons and pions, which allows for both local as 
well as pion-mediated multi-nuclcon interactions, was de- 
veloped in Refs. fl2| in the heavy-baryon limit and the 
main trends for the energy density of symmetric nuclear 
and neutron matter were already reproduced at next-to- 
lcading order in their power counting scheme. 

In applying TFT to nuclear matter in a relativistic set- 
ting, we have to maintain the covariance of the theory 
by using fully relativistic thermal propagators for nu- 
cleons and pions. Keeping the Lorcntz symmetries in- 
tact is of advantage, in particular, for computing the 
scattering and radiation amplitudes with full propaga- 
tors and renormalized vertices. The second goal of our 
work is to maintain self-consistency among the propa- 
gators and the self-energies of the theory, which means 
that the iterations are performed until the Schwinger- 
Dyson equation for the nucleons is fulfilled. If a firm per- 
turbativc expansion (with power-counting rules) exists, 
a self-consistent approach generates higher-order terms 
in every iteration, and therefore is not really necessary. 
Nevertheless, using the leading-order term in the chiral 
expansion of the Lagrangian, we will check the impact 
of these higher-order terms in a self-consistent solution. 
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Third, wc shall address the role of the two- and three- 
body contact terms which are necessary to achieve nu- 
clear saturation in isospin-symmctric nuclear matter and 
we shall discuss some of the key nuclear observables such 
as the symmetry energy and the compressibility. 

This paper is structured as follows. In Sec. HI] we dis- 
cuss the chiral Lagrangian. Section IIIII uses TFT to 
compute the pion contribution to the nucleon self-energy. 
Our numerical method and results for the self-energy as 
well as nuclear matter properties are presented in Sec. lIVI 
Our conclusions are collected in Sec. |V] We use natural 
units h = c = ks = 1. Four- vectors are denoted with 
capital letters, for instance P M = (p°,p). 



II. LAGRANGIANS 

Low-energy nuclear dynamics can be constructed on 
the basis of the pion and nucleon degrees of freedom 
starting from a chiral Lagrangian. The interaction La- 
grangian £„n between nucleons and pions is constructed 
such as to reflect the spontaneous chiral symmetry break- 
ing of strong interactions at low energies. Since the 
interactions of Goldstonc bosons must vanish at zero- 
momentum transfer and in the chiral limit (i.e., the pion 
mass ttt-tt — > 0), a low-energy expansion in powers (the 
so-called chiral dimension) of the ratio of the momentum 
or the pion mass over (4tc times) the pion-decay constant 
can be performed. Consequently, the Lagrangian can be 
written as 



r 
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(1) 



where the superscript labels the order of the chiral di- 
mension. The terms in the expansion ((T|) arc constructed 
by introducing the following SU (2) matrix U in flavor 
space 



U = cxp 



1 + ■ TV 



1 

2/f 



(2) 



where r is vector of Pauli matrices in isospin space, tv is 
the isotriplet of pions, and f v the pion decay constant. 
The leading-order term is given by [13| 



4^ = ^1*7"^ 



(3) 



where tp is the nucleon field, ip = V* 70j m is the nucleon 
mass, and gA is the axial-vector coupling. The physical 
value of gA is determined from neutron beta decay and 
is given by g A = 1.2695 ± 0.0029 
derivative, 



D M is the covariant 



(4) 



where is the so-called chiral connection which couples 
an even number of pions to the nucleon and is defined as 

r l t = i{S\d lt t} = i@d ll z + tdptf) 
i 



■ r ■ (tt x d^Tv) + 



(5) 



with £ = y/U. The La grangian ([3]) also includes the axial- 
vector current which couples an odd number of pions 
to the nucleon 

u„ = i[£},dS = i$d£ - £0 M £t) = ~t • 0„7r + . . . . 

(6) 

Keeping only the lowest-order term in the chiral La- 
grangian ([3]) the pion-nucleon Lagrangian reads @] 
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where we have neglected the Weinberg- Tomozawa contri- 
bution arising from the chiral connection ([5]). The chiral 
onc-pion interaction term in Eq. takes the light dy- 
namical degrees of freedom, i.e., pions, explicitly into 
account. The complete interaction Lagrangian includes, 
in addition, the four-fermion and, in addition, the six- 
fermion contact interaction, respectively, which become 
important at high density. Thus, the Lagrangian of the 
system can be written as 



£ = £ 
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where the first term is the Lagrangian of the non- 
interacting system and the last two terms correspond 
to the two-body and three-body interactions. Figure Q] 
shows the three distinct contributions to the self-energy 
of a nucleon originating from the tree-level vertices in the 
Lagrangian jS]). 



III. PION CONTRIBUTION TO THE 
SELF-ENERGY 

A. Leading-order contribution 

In this section we evaluate the one-pion contribution 
to the nucleon self-energy, corresponding to the left di- 
agram in Fig. [T] within the imaginary-time formalism. 
The free covariant propagator of the nucleons in energy- 
momentum space is given by 



S (K) 



At 



A fe7o 



ik n + a* - Ek ik n + A* + Ek 



(9) 



where the zeroth component of the four-momentum 
= (fc°,k) takes discrete values, fc° = ik n = (2n + 
l)iriT, n £ Z and T is the temperature. are projec- 
tors onto positive (+) and negative (— ) energy states, 



1 



A ^ = 2E~ k [Ek ± ( " ' k + m7 ° )] 



(10) 



where E 2 



2 + m 2 is the dispersion relation for non- 
interacting nucleons, m is their mass, and a = 70 7- The 
free pion propagator is given by 



D (Q) 
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2oj„ 
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(11) 
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FIG. 1. The three diagrams contributing to the self-energy of the nucleon. The first one represents the chiral one-pion exchange 
contribution to the nucleon self-energy. Solid lines refer to nucleons, dashed lines to pions and dots to pion-nucleon vertices. 
The second diagram is the contribution from the two-body contact interaction, with the square vertex representing the two- 
body scattering matrix. The third diagram is the contribution from the three-body interactions with the hexagonal vertex 
representing the three-body scattering matrix. In each diagram the external legs do not contribute to the self-energy. 



where the zeroth component of the four-momentum = 
(q , q) takes discrete values q° = ioj n = 2niriT, n € Z, 



and ui 2 q 



q 2 + m 2 is the dispersion relation for non- 



interacting pions, with m v being the free pion mass. In 
terms of the free propagators ([9]) and pip the one-pion 
exchange contribution to the nucleon self-energy reads 



summation can be arranged according to these under- 
lying processes, but it is more convenient to separate 
the self-energy into the "vacuum" and "thermal" parts 
£(P) = Eo(-P) + ^t(P), where the vacuum part is given 
by 



£(P) 



4/2 



Arp 



(12) 



where P M = (p°,p) = (ip n ,p)- We substitute the 
propagators ^ and (fTT) into Eq. (JT3J) and carry out 
the summation over the fermionic Matsubara frequency 
ik n . In general, this sum generates physically distinct 
processes involving all possible combinations of bosons 
and fermions and their antiparticles. The result of the 



3g 2 A f d 3 k 1 
4/2 J (2tt) 3 2uj pk 
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and the thermal part is given by 



35i 



d 3 k i 

(27t) 3 2oj pk 
1 



1 



/' 



1 



Jpk IPn + V- - Ek + u pk _ 
1 



75 (i 6 - 



where uj 2 ,. 



75 (f 

-k) 2 - 



/i + 

-t) 
■t) 



E k - OJpk iPn + M - 
A fc7Q 

ip n + fl- E k - UJ p k 

A fc7o 



E k 



75 

A fe7o 



^)A+ 7o7 5(i 6 



E k - 



n F (E k - n) 



ip n + H + Ek- OJpk 

A fc7o 



IPn 



E k 



k + LOpk IPn + A* + Ek + OJpk 



75 { 



75 { 



*) 



k =ip„+ij pk 



riB{ojpk) 



(14) 



mt and nmnix) = [exp(x/T)± 1] 1 are the Fermi/Bose distribution functions. The retarded 



self-energy is obtained by analytical continuation, i.e., ip n — > po + i0 + . We have verified that in the case of a Yukawa 
interaction Eq. (|14p transforms to the well-known expression for the self-energy of a fermion in finite-temperature 
quantum field theory 0, [l4[ . 

At sufficiently low temperature the occupation number of anti-particles is so small that we can neglect their 
contribution in Eq. (j 14|) . Furthermore, we assume that there is no macroscopic occupation of pionic modes in nuclear 
matter at any temperature and density of interest; therefore, we also drop the contributions proportional to the 
bosonic occupation numbers. The remaining contribution arises from the pole at ko = Ek — [i and we arrive at 
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where we have defined the quantities s 

1 



3o 

q 



[{p -k y + (p-kY} 



m/2E k )(P — K) 2 and Q M = (qo, q), with components 
1 



E k 



(p - k) • k(p - k ) , 



1 

2E k 



[(P - K) 2 k + 2(p - k) ■ k(p - k)] + (po - ko)(p - k) 



(16) 
(17) 



Later on we will enforce self-consistency in evaluating the self-energy. This requires a Lorcntz decomposition of the 
self-energy, which in the most general case is given by 



E(P) - E S (P) + 75 E ps (P) + 7 M S M (P) + 757 M 5> (P) + ^E^(P) 



(18) 



The requirements of parity conservation, translational and rotational invariance, as well as time-reversal invariance, 
reduce this most general decomposition to the following form 



E(P) = E S (P) + 7o Eo(P) + 7 • PE„(P) • 



(19) 



Equation (j!5p can now be projected onto its Lorcntz components by multiplying it with l,7o, and 7, and taking 
the trace over the 7-matrices. Keeping only the thermal part of the self-energy (and dropping the index T on the 
sclf-encrgies), this leads us to the following decomposition coefficients 



%| r rf 3 k 1 

%s[F) - APJ (2^2^ 



So(P) 



pk 

d 3 k 1 



1 



Po - k - LOpk + i0+ p - fc + Lu pk + i0+ 
1 1 



4/2 J {2vf2uj vk 
l((Po - fco) 2 + (P - k) 2 ) + i-(p - k) • k( Po - ho) 



Po - fc - uj pk + i0+ p a - fc + + i0+ 



n F {ko) — {P-Kf 
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ph ^ F) ~ 4f 2 J (2^2^ 
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n F (k ) 



— [(P - K) 2 k + 2(p - k) • k(p - k)] - (p - fco)(p - k) \- p , 



(20) 



(21) 



(22) 



where p = p/|p| is a unit vector. Equations (f2TJ|) (l22t are our final result for the chiral one-pion-exchange to 
the nucleon self-energy. It is evident that the other terms in Eq. (|14p . which could become important at higher 
temperatures and densities, can be evaluated in a completely analogous way. 
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B. Further approximations 

For numerical computations the factor (s + (ff) in 
Eq. (|T5|) can be simplified. We start by rewriting the 
expression 

2E k Q = 2E k (q j - q ■ 7) = M [(Po - k Q ) 2 + (p - k) 2 ] 7o 
+ 2(P - K) ■ Kf - (P 2 - K 2 )t 
-2/i(p -fc )(p-k)- 7 . (23) 

For the densities and temperatures of interest, nuclcons 
are constrained to the vicinity of their Fermi surface, 
therefore the momentum of the nucleon can be expressed 
as p = p F n + Sp. Here, p^ is the nucleon Fermi momen- 
tum, n is a unit vector, and 5p is the residual momentum, 
with \Sp\ <C Pf- Furthermore, the relativity parameter 
x = pf/jtl is small as well; numerically, we have 



where n is the density of the systems, no = 0.16 fm~ 3 
is the nuclear saturation density. Then, the nucleon en- 
ergy is E p « 771(1 + p 2 /2m 2 ) and the chemical potential 
/1 « 771(1 + p 2 F /2m 2 ). This implies that po = E p — u « 
(p 2 — p 2 p)j2m « xh ■ 5p is small compared to |p| « p^- 
Therefore, we can replace (P — K) 2 ss — (p — k) 2 . 
With these approximations we obtain 

2E k $ « 2m $ w p(p - k) 2 7o - 2(p - k) • kf 
+ (p 2 -k 2 )^-(p 2 -k 2 )(p-k). 7 
«(p-k) 2 [(p + p) 7o -p- 7 ] (25) 

where we used the fact that po — ko ~ (p 2 — p 2 F )/2m — 
(k 2 -p|)/2m= (p 2 -k 2 )/2m. 

Therefore, the factor s + $ in the integrand of the self- 
energy (|15p reads now 



s + 



(P~k) s 
2m 



[?77. + (po + p) 7 - 7 ' P] 



(26) 



x w 0.28 I — 

,"0 



1/3 



(24) 



Separating the contributions from particles and anti- 
particles in the boson propagator, with these approxima- 



tions the self-energy reads 

E(P) w [m + (po + Mho - 7 • P] [Mi*) + CT -( P )] 

= [to + (p + m)7o - 7 ' P] <r{P) , (27) 

where 

3g 2 A f d 3 k 1 (p-k) 2 



8m/ 2 J (2tt) 3 2w pfe Po ~ T + in 



n F (k ) 
(28) 

Therefore, the coefficients of the Lorentz decomposition 
of the self-energy can be expressed in terms of a(p), 

£ s (P)«toct(P), Eo(P)w(po+MM^). 
E„(P) w -a(P) . (29) 

The real part of the self-energy can now be computed 
with the help of the Dirac identity. We obtain by explic- 
itly evaluating the Cauchy principal value of Eq. (l28l) 



Re[a(P)] = Re[(i + (P) + a-(P)] 



32tott 2 / 2 



iiki k 2 " F(fco) 

11 2|p||k| 



-1 x-x 
(30) 
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FIG. 2. (Color online) Real part of the nucleon self-energy 
as a function of momentum p at n = 0.1 fm -3 and T = 0. 
The dotted (blue) line is the lowest-order result and simul- 
taneously the starting point of the iteration. The full (red) 
line represents the final self-consistent result. The other lines 
correspond to 2,3,4,5 (from bottom to top) iteration steps. 



where x = [p 2 + k 2 + m 2 - (p - /c ) 2 ]/(2|p||k|). 



C. Enforcing self-consistency 

The self-consistency of the numerical computation of 
the self-energy is achieved by replacing the free nucleon 
propagator in Eq. (fT2j) by the full nucleon propagator, 
determined by the Schwingcr-Dyson equation 

S-\P) = Sv\P) - S(P) - [(po + Mho - 7 • P] 

x [1 - Reo-(P)] - m[l + Rccr(P)] + ImS(P), 

(31) 

where the free nucleon propagator is given by Sq 1 (P) = 
(po + m)7o — 7 • p — ?ti and we have used Eq. (j27|) for the 
nucleon self-energy. The roots of det5 _1 (P) determine 
the excitation spectrum of the system, 

pi = E* p - /x* , (32) 

where 

E* = ^p* 2 +m* 2 , (33) 

m* =m[l + Rea(P*)} , (34) 

p* = p[l - Re<7(P*)] , (35) 

M* = m[1 - Re«r(P*)] . (36) 

We achieve self-consistency for the self-energy by replac- 
ing the free quantities m, p, k, and fi in the integrand 
(but not in the integration measure) of Eq. (f3T)|) by the 
corresponding rcnormalized m* , p* , k* , and fi* . In prac- 
tice, we start by computing (|30j) with the free quanti- 
ties, which define the renormalized quantities (|33[) - (|36[) 
to first order in the iteration process. We repeat the 
previous step until convergence is reached. 



■-■ T = 20MeV 
T = 0MeV 




P [fin"'] 

FIG. 3. (Color online) Dependence of the one-pion exchange 
contribution to the nucleon self-energy in isospin-symmetric 
nuclear matter on the momentum at density n = 0.1 fm~' ! . 
The solid line shows the real part of the self-energy for tem- 
perature T = and the dashed-dotted line for T = 20 MeV. 



IV. RESULTS 

The energy per nucleon, excluding its rest mass to, is 
defined as 

^ = T + W-m, (37) 

where T is the kinetic energy, U the single-particle poten- 
tial and N is the particle number. In terms of the renor- 





BHF SCGF Viow-k ChPT 


to [MeVfm 3 ] 
ts [MeV fm 3+3a ] 


-153 -311 -438.1 -836.052 
2720 3670 6248 8048.29 



TABLE I. Parameters to and t^ defining the contact inter- 
action of Eq. (|40p as obtained from a fit to the saturation 
point n = 0.16 mi -3 and E/N = -16.0 MeV at a = 0.5 for 
different models of nuclear matter [l5| . 



malized quantities ([3"4" |) - (j31))) the kinetic energy reads 

IP 



T 



d\ P \p 2 



E* 



„ \p\)n P (E;-n*), 
p 



(38) 

where g T is the isospin degeneracy factor (g T = 2 in 
isospin-symmctric nuclear matter and g T = 1 in neutron 
matter). The single-particle potential, IA, is given by 



U 



fJr 



d\p\p 2 



2n 2 n 

x a(P)n F (E; - (i*) 



m 

—m + E, 

E* 



p 



E* 



(39) 



As we consider only spin-unpolarized matter, the spin 
summation has been carried out in Eqs. ([55)1 and 



A. Symmetric nuclear matter 

The Fock diagram evaluated in the previous section 
represents chiral onc-pion exchange between two nucle- 
ons, therefore its contribution to the potential energy is 
attractive. In order to achieve saturation in symmetric 
nuclear matter, we need to add a potential which includes 
the self-energy contribution arising from the two- and 
three-nucleon contact interactions, shown in Fig. [TJ This 
contribution can be expressed in the form of a Skyrme- 
likc interaction, given by 



V c (p) = lt n+^t 3 n 1+a , 



(40) 



where n is the density and to, £3, and a arc parameters. 
For a fixed value of a (typically a = 0.5), we fit to and 
tj, in such a way that Eq. (|37p yields the empirical satu- 
ration point for symmetric nuclear matter. For a = 0.5, 
the fit results in an attractive two-body contact interac- 
tion (t < 0) and a repulsive three-body one (t 3 > 0). 
The numerical values for to and ts are given in Table fl] 
We note that the value of a controls the stiffness of the 
equation of state at high densities, with larger values cor- 
responding to a stiffer equation of state. The particular 
value a = 0.5 is consistent with the density-dependent 
part of the Skyrme interaction. 

Once the parameters are fixed, we compute the nuclear 
compressibility modulus, defined as 



K = 9n 



,d 2 (E/N) 



dn 2 



(41) 




FIG. 4. (Color online) Dependence of the energy per nucleon 
in symmetric nuclear matter on density at T = (solid red 
line). The chiral one-pion exchange contribution to the po- 
tential energy is shown by the dashed-dotted (green) line and 
the contribution from contact terms by the dashed (blue) line. 
For comparison, we show the one-pion exchange contribution 
(dots) as well as the full result (squares) from calculations 
using a low-momentum potential Viow-k. 



We find K sa 267 MeV, which is in agreement with the 
empirically known value (K « 250 ± 25 MeV). 

The self-energy due to chiral onc-pion exchange in 
isospin-symmctric nuclear matter is shown in Fig. [5] at 
T = and n = 0.1 fm~ 3 . As one observes, the self- 
consistency procedure converges rapidly. In order to 
achieve a relative accuracy < 10~ 6 one needs about 10 
iterations, the exact number of iteration depending on 
the density. For small densities the number of required 
iterations is small, while for large densities the number of 
iterations needed to achieve convergence is larger. While 
the self-energy does not tend to zero for large external 
momenta, the result for the energy of the system is con- 
vergent due to an additional integration over a Fermi 
distribution function. 

In Fig. [3] we show the self-energy at zero temperature 
in comparison to its form at T = 20 MeV. The change 
with temperature can be seen to be rather moderate. In 
the non-zero temperature case the self-energy is larger for 
low momenta than at zero temperature, while for higher 
momenta the variation of temperature does not affect 
the result very much. This is due to the fact that we 
deal with comparatively low temperatures, therefore the 
temperature is a relevant scale only at low momenta. 

It can be seen from the expression of the self-energy in 
Eq. (f2~T)) and from the definition of the self-energy coef- 
ficients (|29|). that the vector self-energy is substantially 
smaller than the other contributions to the self-energy, 
as expected. Our numerical result for the real part of the 
reduced self-energy is in good agreement with the one 
quoted by Ref. [|| in their zero-temperature calculations. 

The energy per nucleon in isospin-symmctric nuclear 
matter at T = is shown in Fig. @] along with its 
potential-energy components which are split into the con- 
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tribution of onc-pion exchange and the contribution from 
two- and three-body contact interactions. The one-pion 
exchange contribution is linear in the density, as one 
would expect from a mean-field type contribution and 
is always negative, as it reflects the attractive interaction 
between nucleons. The two-body contact interaction is 
attractive, while the three-body interaction is repulsive. 
As stated above, the parameters of these interaction were 
chosen to fit the saturation point of nuclear matter. It is 
remarkable that quantitatively the contact interactions 
make a substantial contribution at all densities of inter- 
est, therefore one may conclude that the relativistic one- 
pion exchange alone cannot account for the attraction in 
isospin-symmetric nuclear matter and cannot lead to the 
correct saturation properties in combination with only 
repulsive three-body forces acting at high densities. In 
the limit of vanishing density the kinetic energy becomes 
dominant and we recover the free Fermi gas behavior of 
the energy. 

Figure 2] also shows a comparison of the energies per 
nucleon in isospin-symmetric nuclear matter in our case 
and in the model which uses the Viow-k interaction. The 
latter one is a soft, phase-shift equivalent interaction 
which applies below some momentum cut-off. As can 
be seen, this model also requires a potential-energy con- 
tribution from two- and three-body contact interactions 
in order to achieve saturation of nuclear matter. The 
models agree well (per fit of the contact interactions) for 
the energy per nucleon, but in both cases a substantial 
portion of the attractive energy has to be included in 
the phenomenological two-body contact term, in order 
to reproduce the phenomenological features of isospin- 
symmetric nuclear matter. 

The compressibility in our model is K = 267 McV, 
which is in good agreement with the one predicted by 
the VW_k interaction (K = 258 McV). Table fl] com- 
pares the parameters of the two- and three-body contact 
interactions required to fit the nuclear saturation point in 
isospin-symmetric nuclear matter for several models. In- 
deed, all non-relativistic models overbind nuclear matter 
at the saturation density, therefore repulsive three-body 
forces are mandatory to cure this problem. It is seen that 
the values used in the literature and those in our fits are 
more or less of the same order of magnitude and, for in- 
stance, the ratios between the parameters in our model 
to those of the Mow-k interaction are ~ 2 for to and 4/3 
for t 3 . 



B. Pure neutron matter 

The self-energy of neutrons in pure neutron matter is 
shown in Fig. [5] Compared to isospin-symmetric nuclear 
matter, the self-energy in neutron matter is smaller. This 
is due to the fact that in pure neutron matter the one- 
pion exchange involves only the 7r°-meson. At non-zero 
temperature, as in the case of isospin-symmetric nuclear 
matter, we observe that, for low momenta, the contri- 
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FIG. 5. (Color online) Dependence of the one-pion exchange 
contribution to the nucleon self-energy in pure matter on the 
momentum, at density n = 0.1 fm -3 . The solid line shows 
the real part of the self-energy for temperature T — while 
the dot-dashed line - for T = 20 MeV. 



bution of the self-energy at higher temperature is larger, 
but the difference among the two cases tends to disappear 
with increasing momentum. 

The energy per particle in pure neutron matter in 
shown in Fig. [5] Now we fit the parameters at two val- 
ues of the density, namely at saturation and at twice 
the saturation de nsity , to the quantum Monte-Carlo cal- 
culations of Ref. [lfj, which are based on the Argonne 
AV8' two-body and Urbana IX three-body forces. The 
attractive contribution from chiral one-pion exchange is 
smaller in this case, due to the smallness of the neutron 
self-energy. Again, we see that the contribution from con- 
tact interactions is substantial in order to get an equa- 
tion of state compatible with microscopic calculations. 
The fit parameters in the present case are to = —346.64 
and £3 = 5255.08 and they differ from those in isospin- 
symmetric nuclear matter due to the isospin dependence 
of the nuclear force. 

Finally, wc report our result on the dependence of 
the symmetry energy S on density in nuclear matter in 
Fig. [7J The empirical value of S at the saturation point 
is S(no) = 32.5 McV and we find the value S(n ) = 35.1, 
which is a consequence of the fit to the neutron matter 
equation of state. 



V. CONCLUSIONS 

In this work we have combined the methods of TFT 
and chiral Lagrangians to compute the self-energy of nu- 
cleons to leading order in the chiral expansion. In doing 
so, we have maintained the covariancc of the pion and nu- 
cleon propagators and we have imposed self-consistency 
by solving a Schwingcr-Dyson equation for the nucleon 
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FIG. 6. (Color online) Dependence of the energy per particle 
in pure neutron matter on density at T = (solid red line). 
The contribution from chiral one-pion exchange is shown by 
the dash-dotted (green) line and that from the contact inter- 
actions by the dashed (blue) line. 
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FIG. 7. Dependence of the symmetry energy of nuclear matter 
on density at T — 0. 

self-energy. Medium- and short-range interactions be- 
tween nucleons are represented by two- and three-body 
contact interactions which are chosen to have the struc- 
ture of the corresponding terms in Skyrme density func- 
tionals. These were then fitted to reproduce the nuclear 



saturation in the case of isospin-symmetric nuclear mat- 
ter. 

We find that the relativistic one-pion exchange con- 
tributes about half of the saturation energy for isospin- 
symmetric nuclear matter, the remainder being provided 
by two- and three-body contact interactions. By fit- 
ting the contact interactions to the saturation point for 
isospin-symmetric nuclear matter, our model produces 
reasonable values of the compressibility modulus. A com- 
parison between our results and those obtained from a 
Mow-k interaction shows a good agreement for the attrac- 
tive contribution to the energy per nucleon for isospin- 
symmetric nuclear matter. 

Our approach has also been applied to pure neutron 
matter with similar results: the relativistic one-pion con- 
tribution needs to be supplemented with two and three- 
body contact interactions in order to match an equation 
of state derived from microscopic calculations. 

Clearly, to obtain a consistent phenomenology of both 
isospin-symmetric nuclear and pure neutron matter, one 
needs to introduce contact interactions which account for 
the short-range attraction at the two-body level, as well 
as for three-body repulsion. In this respect further steps 
might be undertaken to resolve the relativistic dynam- 
ics of pions by including higher-order terms in the chiral 
expansion and incorporating A-isobar excitations. 

Methodologically, our approach differs from similar 
works since we address the nuclear equation of state with 
a chiral Lagrangian keeping a relativistic framework, and 
at the same time we impose self-consistency by solving a 
Schwingcr-Dyson equation. 

The in-medium electromagnetic and weak interactions 
of nucleons can be computed in a relativistically covari- 
ant manner starting from the self-consistent propagators 
derived above. 
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